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ABSTRACT 

Histone deacetylases (HDACs) and lysine acetyl- 
transferases (KATs) catalyze dynamic histone 
acetylation at regulatory and coding regions of 
transcribed genes. Highly phosphorylated HDAC2 
is recruited within corepressor complexes to 
regulatory regions, while the nonphosphorylated 
form is associated with the gene body. In this 
study, we characterized the nonphosphorylated 
HDAC2 complexes recruited to the transcribed 
gene body and explored the function of HDAC- 
complex-mediated dynamic histone acetylation. 
HDAC1 and 2 were coimmunoprecipitated with 
several splicing factors, including serine/arginine- 
rich splicing factor 1 (SRSF1) which has roles in 
alternative splicing. The co-chromatin immunopre- 
cipitation of HDAC1/2 and SRSF1 to the gene 
body was RNA-dependent. Inhibition of HDAC 
activity and knockdown of HDAC1, HDAC2 or 
SRSF1 showed that these proteins were involved 
in alternative splicing of MCL1. HDAC1/2 and 
KAT2B were associated with nascent pre-mRNA 
in general and with MCL1 pre-mRNA specifically. 
Inhibition of HDAC activity increased the occu- 
pancy of KAT2B and acetylation of H3 and H4 of 
the H3K4 methylated alternative MCL1 exon 2 nu- 
cleosome. Thus, nonphosphorylated HDAC1/2 is 
recruited to pre-mRNA by splicing factors to act 
at the RNA level with KAT2B and other KATs 
to catalyze dynamic histone acetylation of the 
MCL1 alternative exon and alter the splicing of 
MCL1 pre-mRNA. 



INTRODUCTION 

Lysine acetyltransferases (KATs) and histone deacetylases 
(HDACs) catalyze dynamic acetylation of proteins, 
including histones, associated with transcribed DNA (1-7). 
KATs often have transcriptional coactivator activity, 
increasing the level of acetylated histones and enhancing 
transcription when recruited to a gene promoter by a tran- 
scription factor (1,8-10). HDAC1 and HDAC2 are present 
in multiprotein corepressor complexes (Sin3, NuRD, 
CoREST), which are recruited to regulatory regions by tran- 
scription factors (11,12). Phosphorylation of HDAC1 (at 
S393/421/S423) and HDAC2 (at S394/422/S424) is 
required for the formation of these corepressor complexes, 
which are recruited to the regulatory regions of transcribed 
genes (11-14). In contrast, the nonphosphorylated HDAC2 
is associated with the body of transcribed genes (12). 
Although the unmodified HDAC2 is more abundant 
than highly phosphorylated HDAC2, it is the highly 
phosphorylated form that is preferentially crosslinked to 
chromatin with formaldehyde or cisplatin (11). However, 
through the use of a dual crosslinking chromatin 
immunoprecipitation (ChIP) assay, all isofomis of 
HDAC1 and HDAC2 could be mapped along regulatory 
and coding regions of transcribed genes, with the unmodi- 
fied HDAC2 being associated with the coding region (12,15). 
Yet, it remains to be determined which proteins the unmodi- 
fied (nonphosphorylated) HDAC1 and HDAC2 interact 
with when targeted to the body of transcribed genes. 

Recent studies suggest a role for HDAC1 and HDAC2 
in alternative splicing (16,17). Approximately 95% of 
human multiexon genes generate alternatively spliced 
transcripts, giving rise to mature mRNA isoforms 
coding for functionally different proteins. Most of these 
splicing events are regulated in a tissue- and/or devel- 
opmental stage-specific manner or in response to naturally 
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occurring external stimuli (18-20). Pre-mRNA splicing is 
a cotranscriptional process, which is regulated by RNA 
polymerase II (RNAPII) elongation rate, as inclusion of 
an alternative exon occurs only if splicing components 
have time to interact with the nascent RNA before its 
3'-end cleavage and release (21-25). 

There is emerging evidence that histone modifications 
and chromatin structure influence splicing, and vice- 
versa pre-mRNA splicing itself influences chromatin 
organization (23-30). Several studies have suggested or 
demonstrated a correlation between local increase in 
histone acetylation and exon skipping; however, the mech- 
anisms involved are poorly understood (16,17,29,31). The 
skipping of NCAM exon 18 upon membrane depolariza- 
tion of neuronal cells was linked to localized increased 
histone H3 acetylation (H3K9ac) and could be replicated 
by the HDAC inhibitor trichostatin A (TSA) (29). In HeLa 
cells treated with the pan-HDAC inhibitor, sodium 
butyrate, splicing-sensitive exon-arrays detected a change 
in the splicing pattern of ~700 genes (16). The levels of 
serine/arginine-rich (SR) proteins and other proteins 
involved in splicing, as well as the acetylation status of 
splicing factors were not altered by 1 5 h of HDAC inhib- 
ition. In the case of the fibronectin (FN1) alternative exon 
25 (also named EDB), histone H4 acetylation increased 
rapidly following HDAC inhibition and reached its 
maximal level after 6-9 h. An increased exclusion of exon 
25 could be detected within this time frame, if cells were 
treated with 5,6-dichloro-l-p-D-ribobenzimidazole (DRB), 
a reversible inhibitor of RNAPII, such that analysis would 
be limited to de novo synthesized pre-mRNA (16). It was 
shown that HDAC inhibition decreased the association of 
one of the SR proteins, SRSF5 (also known as SRp40) with 
the FN1 gene, including but not restricted to exon 25. 
siRNA-mediated knockdown of HDAC1, but not 
HDAC2, resulted in exon 25 skipping, suggesting that 
HDAC1 is primarily involved in the splicing regulation of 
this gene (16), but the mode of action of HDAC1 remains 
unclear. In neuronal cells, it was suggested that HDAC2 
association with the splicing regulator Hu proteins 
enhanced the localized histone acetylation at the alternative 
exons of NF1 and FAS genes, an event that was correlated 
with a localized increased elongation rate and the exclu- 
sion of these exons. In vitro, the HuR proteins inhibited 
HDAC2 activity (17). It was proposed that HuR proteins, 
cotranscriptionally recruited to their target RNA se- 
quences, inhibit HDAC2 activity through a 'reach back' 
interaction with chromatin (17). However, the association 
of HDAC2 with the alternative exons of NF1 and FAS 
genes was not directly demonstrated. Thus, a correlation 
between histone acetylation and skipping of alternative 
exons was observed, but the recruitment and distribution 
along the body of the transcribed genes of HDACs 
and KATs, the enzymes catalyzing histone acetylation, 
were not addressed in these studies. 

In this study, we demonstrate that the serine/arginine- 
rich splicing factor 1 (SRSF1, also known as ASF/SF2) 
coimmunoprecipitated with HDAC1 and nonpho- 
sphorylated HDAC2. Since this is a SR protein with a 
major role in alternative splicing (32), we explored the 
role of HDAC1 and HDAC2 in splicing using the 



human myeloid cell leukemia sequence 1 (MCL1) gene 
as a model gene. The MCL1 gene undergoes alternative 
splicing of exon 2 and produces a protein that either 
prevents or supports cell death depending on the splice 
isoform. The long form MCL1L is an anti-apoptotic 
protein, while the short form MCL1S is pro-apoptotic 
(33). We have analyzed the effects of HDAC pan- and 
Class I specific inhibitors on the splicing of the MCL1 
gene. We show that HDAC inhibitors enhanced histone 
acetylation over exon 2, an event that paralleled exon 2 
exclusion. We also show, for the first time, that HDAC1 
and HDAC2 are recruited to the MCL1 pre-mRNA to 
catalyze, in concert with KATs, dynamic histone acetyl- 
ation of the exon 2 nucleosome. 

MATERIALS AND METHODS 

Cell cultures and treatments 

HEK293, Flp-In 293 stable cells expressing HDAC2-WT - 
V5 or HDAC2-3S/A-V5 and MCF7 cells were cultured in 
Dulbecco's modified Eagle's media (Gibco), and HCT116 
cells were cultured in McCoy's 5A media (Sigma), supple- 
mented with 10% fetal bovine serum (FBS), 100 units/ml 
penicillin, lOOug/ml streptomycin, and 250ng/ml ampho- 
tericin B, and were maintained at 37°C in a humidified 
atmosphere containing 5% C0 2 . Cells were treated with 
HDAC inhibitors, TSA (250 nM) (Sigma), apicidin 
(150 nM) (Sigma) or butyrate (2mM) (Sigma) for 
indicated time periods. When indicated, subconfluent 
cells (~80-90%) were serum-starved for 48 h and then 
either left untreated or treated with TPA (12-O-tetradeca- 
noylphorbol-13-acetate) (100 nM, Sigma) for 30 or 60min. 
In inhibition studies, serum starved cells were pretreated 
with either TSA or apicidin for 30 min followed by treat- 
ment with phorbol ester, TPA (0, 30 or 60 min). 

Mass spectrometry 

HDAC1 or HDAC2 complexes were immunoprecipitated 
from HEK293 cells or Flp In 293 stable cells expressing 
WT-HDAC2-V5 or HDAC2-3S/A-V5, with anti- 
HDAC1, anti-HDAC2 (Affinity BioReagents) or anti-V5 
(Abeam) antibodies. When mentioned, nuclear extracts 
from Flp In 293 stable cells expressing WT-HDAC2-V5 
or HDAC2-3S/A-V5 were treated with RNase A (5 ug of 
RNase A for 20 min at 37°C) prior to immunoprecipitation 
with anti-V5 antibodies. Covalent immobilization of the 
antibodies onto the surface of Dynabeads® Protein G 
(Invitrogen) was used for coimmunoprecipitation of 
intact protein complexes. The immunoprecipitated frac- 
tions were eluted from the Dynabeads with 1% SDS/ 
0.1 M NaHC0 3 . The eluted fractions were vacuum dried 
and washed with 100 mM NH4HCO3 and iodoacetamide. 
After lyophilization, the fractions were digested with 
trypsin for 16 h at 37°C. The nano-liquid chromatography 
and tandem mass spectrometry were performed on the 
trypsin-digested samples as described previously (34). The 
MSDB, version 20060831, database was searched using the 
Global Proteome Machine (http://www.thegpm.org) 
search engine to identify the peptide sequences. 
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Immunoprecipitation and immunoblotting 

Cells were harvested and lysed in cold lysis buffer (50 mM 
Tris-HCl, pH 8.0, 150mM NaCl, 1 mM EDTA, 0.5% 
NP-40) containing phosphatase and protease inhibitors 
(Roche), and immunoprecipitations were done as 
described earlier (11). In brief, 500 ug total cellular or 
nuclear extracts were incubated with 3.0 ug of anti- 
HDAC1 (Affinity Bioreagents), anti-HDAC2 (Affinity 
Bioreagents), anti-SRSFl (Santa Cruz Biotechnology), 
anti-acetyl lysine (Cell Signaling Technology), anti-H3 
(Millipore) or anti-RNAPIIS2ph (Abeam) antibodies 
overnight at 4°C. An amount of 40 ul of protein A/G 
UltraLink resin (Pierce) were added and incubated for 
3h at 4°C. The beads were then washed four times with 
ice-cold lysis buffer. Immunoprecipitation with isotype 
specific nonrelated IgG was performed as negative 
control. One-third of the fractions immunoprecipitated 
by anti-HDACl or anti-HDAC2 antibodies and one-half 
of the fractions immunoprecipitated by anti-RNAPIIS2ph 
or anti-SRSFl antibodies were analyzed by immunoblot- 
ting, while equivalent volumes of lysate (Input) and 
immunodepleted fractions were analyzed, corresponding 
to 20 (for anti-HDACl or anti-HDAC2 antibodies) or 
25 ug (for anti-RNAPIIS2ph or anti-SRSFl antibodies) 
of lysate proteins. When mentioned, cellular extracts 
were treated with RNase A (400ug/ml of RNase A for 
30 min at 37°C) before immunoprecipitation reaction. 
Immunochemical staining was performed with rabbit 
polyclonal antibodies against human HDAC1 (Affinity 
BioReagents), HDAC2 (Affinity BioReagents), HDAC2 
(phospho S394) (Abeam), SIN3A (Affinity BioReagents), 
RCOR1 (Abeam) or mouse monoclonal antibodies 
against HDAC2 (Millipore) or RBPP4 (Abeam). 

HDAC activity assay 

HDAC activity assay was performed with the Fluor-de- 
Lys® HDAC fluorometric activity assay kit (Enzo life 
sciences) following the manufacturer's instructions. 
SRSF1 complex was immunoprecipitated from 500 ug of 
HCT116 cell lysates with 5ug of mouse monoclonal anti- 
SRSFl antibody (Santa Cruz Biotechnology). As a 
negative control, immunoprecipitation with isotype-spe- 
cific nonrelated IgG was performed. An amount of 40 ul 
of protein A/G UltraLink resin (Pierce) were added and 
incubated for 3h at 4°C. The beads were washed three 
times with IP buffer (50 mM Tris-Cl pH 8.0, 150mM 
NaCl, ImM EDTA, 0.5% NP-40) and twice with the 
HDAC activity assay buffer before used for the HDAC 
activity assay. For the assay, the beads treated or not with 
1.0 uM TSA were incubated with 150 uM Fluor de Lys® 
Substrate for 30 min at room temperature with rocking. 
After that, the developer I solution containing 1.0 uM 
TSA was added and aliquots were incubated for another 
30 min to stop the reactions. For boiling control, after 
immunoprecipitation, the beads were boiled for 15 min 
and the assay was performed as mentioned above. The 
fluorescence signal was measured using fluorometric 
plate reader (Spectra MAX GEMINI XS, Molecular 
devices). 



ChIP and re-ChIP assays 

ChIP and re-ChIP experiments were done as previously 
described with an additional protein-protein crosslinking 
step with DSP (dithiobis[succinimidylpropionate]) 
(Thermo Fisher Scientific) (35). Cells were incubated 
with 1 .0 mM DSP for 30 min at room temperature accord- 
ing to manufacturer's instruction (Pierce), followed by 
crosslinking with formaldehyde for 10 min. Dual 
crosslinked chromatin was processed to mono- 
nucleosomes by micrococcal nuclease (MNase) 
(Worthington Biochemical Corporation) digestion (2.5 U 
of MNase/A 2 6o of nuclear suspension) for 40 min at 37°C 
(35). ChlPs were done with anti-RNAPIIS2ph (Abeam), 
anti-HDAC2 (Affinity Bioreagents), anti-SRSFl (Santa 
Cruz Biotechnology), anti-H3K14ac (Abeam), anti- 
H3K9ac (Abeam), anti-acH3 (Millipore), anti-H4K5ac 
(Millipore), anti-H4K8ac (Millipore), anti-acH4 
(Millipore), anti-H3 (Millipore), anti-KAT2B (Abeam) 
or anti-KAT7 (Santa Cruz Biotechnology) antibodies. 
Negative control included performing ChlP/reChIP 
assays with an isotype matched nonrelated IgG. For 
RNase A-treated extracts, DSP and formaldehyde 
crosslinked chromatin was processed to 
mononucleosomes and treated with 400ug/ml of RNase 
A for 30 min at 37°C. The further processing of chromatin 
fragments were performed as previously described (35). 
Input and ChlP/reChIP DNAs were quantified using 
PicoGreen assay. Equal amounts of input and 
immunoprecipitated DNA (1.0 ng) or re-ChIP (0.5 ng) 
DNA were used to perform SYBR Green real-time PCR 
on iCycler IQ5 (BioRad). Primers are described in 
Supplementary Table SI. Enrichment values, calculated 
as previously described (35), are relative to input DNA 
and are the mean of three independent experiments. The 
error bars indicate standard deviation (« = 3). 

ChIP DNA library, sequencing and ChlPSeq 
data analysis 

ChlPSeq libraries for H3K4me3 antibody (Abeam) were 
prepared according to the 5500 SOLiD fragment library 
protocol (Life technologies). Briefly, a quantity of 18 ng of 
ChIP DNA was end-repaired and size selected (100— 
250 bp), followed by dA tailing, ligation of SOLiD 
barcodes and library PCR amplification. Libraries were 
subjected to emulsion PCR (ePCR) and loaded on the 
flowchip for sequencing. Total 18 369 095 and 30049407 
50-base SOLiD sequence reads were generated from 
H3K4me3 ChIP sample and input sample, respectively. 
The sequence reads were mapped on human genome ref- 
erence hgl9 using LifeScope Genomic Analysis Software 
v2.5.1 (Life Technologies) with 2-mismatch setting and the 
best score mapping when multiple mapping occurs. The 
mapping quality was inspected using the Bamstats 
module. More than 83% of reads were mapped. The 
H3K4me3 marker peaks were called using MACS 
package (36) and annotated using CEAS program (37). 
The H3K4me3 signals along the MCL1 gene were 
viewed by the Partek Genomic Suite (Partek Inc. St. 
Louis, USA). 
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siRNA-mediated transient knockdown 

HCT116 cells were transfected at ~30-40% confluency 
with 75 nM of scramble (Nontargeting Pool), human 
HDAC1, HDAC2 or SF2 (SRSF1) ON-TARGET plus 
SMARTpool siRNAs (Thermo Scientific-Dharmacon), 
using Polyplus Interferin siRNA transfection reagent 
(VWR) according to manufacturer's protocol. Forty- 
eight hours after transfections, cells were harvested and 
the knockdown efficiencies were analyzed by immunoblot- 
ting and changes in MCL1 splicing were analyzed by 
radiolabeled PCR ( 32 P). For ChIP assays, cells were dual 
crosslinked after the transfection period and processed as 
described above. 

RNA isolation, reverse transcription PCR and 
radiolabeled- PCR ( 32 P) 

Total RNA was isolated from untreated and treated 
HCT116 cells using RNeasy Plus Mini Kit (Qiagen) ac- 
cording to manufacturer's instructions. Total RNA 
(400 ng) was used as template for synthesis of cDNA 
using M-MLV reverse transcriptase and Oligo dT 
primers (Invitrogen). Radiolabeled 32 P PCR was per- 
formed in a final volume of 15ul containing 7.5 ng of 
cDNA, in the presence of lx PCR buffer (Invitrogen), 
1.5mM MgCl 2 , 0.2mM dNTPs, 0.2 uM of each primer, 
1 unit of Platinum Taq DNA polymerase and 10 nM a 32 P 
dCTP (lOmCi/ml). MCL1 PCR conditions comprised a 5- 
min incubation at 95°C, followed by 35 cycles of (95°C for 
30 s, 56°C for 45 s, 72°C for 45 s) and a final incubation at 
72°C for lOmin. Primers are described in Supplementary 
Table SI. Following PCR amplification, 2.5 jal of 
radiolabeled PCR products were denatured in 80% 
formamide buffer (containing ImM EDTA pH 7.5, 
0.1% xylene xyanol, 0.1% bromophenol blue) and run 
on a denaturing polyacrylamide gel (6%) at 40 W 
(constant wattage). After electrophoresis, the gels were 
dried and exposed to a Molecular Imager TM FX 
(BioRad, Hercules, and CA) and PCR signals were 
quantified. Percentages of the short isoform (MCL1S) 
among total transcripts are presented as percentages of 
exon 2 exclusion. The student's /-test for paired samples 
was used for calculation of statistical significance. 

RNA-CLIP assay (crosslinking and immunoprecipitation 
of RNA-protein complexes) 

Cycling HCT116 cells, treated or not with HDAC inhibi- 
tors TSA (250 nM) or apicidin (150nM) for 2h, were 
either incubated with 1.0 mM DSP for 30min or not, 
followed by irradiation under ultraviolet (UV) light 
(400mJ/cm ) (Stratalinker). Cells were harvested and 
lysed in ice-cold lysis buffer [20 mM Tris-HCl at pH 7.5, 
100 mM NaCl, 10 mM MgCl 2 , 0.5% NP-40, 0.5% Triton 
X100, 0.1% SDS, protease and phosphatase inhibitor 
cocktail (Roche), 80U/ml RNasin (Promega)] with sonic- 
ation. The cellular extracts were treated with a dilute 
cocktail of RNase A/Tl (Ambion) as previously described 
(32). Cellular extracts were precleared with protein A/G 
UltraLink resin (Pierce) and tRNA at lOOug/ml for 2h 
at 4°C. The lysate was incubated with anti-HDACl, 



anti-HDAC2, anti-SRSFl or anti-KAT2B antibodies or 
normal rabbit or mouse IgG on a rotator overnight at 
4°C, followed by addition of protein A/G UltraLink 
resin (Pierce) for 3h at 4°C. Beads were then washed six 
times using the lysis buffer. The beads were resuspended in 
200 ul of lysis buffer and then treated with proteinase K 
(2mg/ml) for 1 h at 55°C. Immunoprecipitated RNA was 
then extracted using Trizol (Invitrogen) and precipitated 
with ethanol. RNA precipitates were resuspended in 20 ul 
of ultrapure water, treated with RQ DNase (Promega) for 
1 h at 37°C. Equal amounts of input and immunopre- 
cipitated RNA were used for cDNA synthesis and 
reversed transcribed using random oligo dT and random 
hexamers with Superscript III reverse transcriptase 
(Invitrogen) following the manufacturer's specifications. 
The resulting cDNA was analyzed by real-time PCR 
with the indicated primer sets in Supplementary Table 
SI. Results are representative of three independently per- 
formed experiments. 

RESULTS 

HDAC1 and HDAC2 interact with splicing factors in a 
RNA-independent manner 

To identify the proteins associated with HDAC1 and 
HDAC2, we established Flp-In 293 stable cell lines ex- 
pressing V5-tagged versions of wild type HDAC2 (WT) 
or a triple mutant HDAC2 (3S/A), which has the three 
CK2 phosphorylation sites mutated from serine to alanine 
(S/A mutation at S394, 422, 424), and immunopre- 
cipitated the HDAC complexes using antibodies 
against the V5 tag (38). We also isolated native 
HDAC 1/2 complexes from HEK293 cells. Proteins 
coimmunoprecipitated with the HDACs were identified 
by mass spectrometry. Table 1 shows that HDAC1 and 
HDAC2 were bound to the Sin3, NuRD and CoREST 
complexes as well as to numerous RNA splicing 
proteins. Nonphosphorylated HDAC2, however, inter- 
acted with RNA splicing factors, but not with components 
of corepressor complexes, confirming previous results that 
HDAC2 phosphorylation is required for corepressor 
complex formation (11-13). Since HDAC1 and HDAC2 
were associated with splicing factors, we determined 
whether these interactions were mediated by RNA. 
Nuclear extracts from Flp-In 293 expressing HDAC2- 
WT-V5 or HDAC2-3S/A-V5 were treated with RNase A 
prior to immunoprecipitation with anti-V5 antibodies. 
The RNase A treatment did not affect the association of 
HDAC2 with RNA splicing factors, showing that the 
interactions between these proteins were RNA-independ- 
ent (Supplementary Table S2). 

Among the RNA-binding proteins associated with 
HDAC1 and HDAC2, we focused our studies on 
SRSF1, which was reproducibly detected in the mass spec- 
trometry analyses. SRSF1 is a sequence-specific RNA- 
binding factor that promotes spliceosome formation by 
binding to exonic splicing enhancers during pre-mRNA 
splicing (39). Interactions between HDAC1 or HDAC2 
and SRSF1 were validated in reciprocal coimmuno- 
precipitation assays using HCT116 cell lysates 
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Table 1. Proteins associated with exogenous wild-type or mutated (3S/A) HDAC2, and with endogenous HDAC1 and HDAC2 
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Except for HNRNPQ, which is coded for by SYNCRIP gene, immunoprecipitated proteins are identified by the genes encoding them. 



(Figure 1A). SRSF1 coimmunoprecipitated with HDAC2, 
but not with either the serine 394 phosphorylated form of 
HDAC2 (HDAC2-S394ph) or the SIN3A, RBPP4 or 
RCOR1 (also known as CoREST), components of 
corepressor complexes (Figure 1A). Thus, the data 
suggest that HDAC2 phosphorylation specifies its incorp- 
oration into corepressor complexes at upstream promoter 
elements of genes, while lack of phosphorylation directs 
HDAC2 association with complexes involved in pre- 
mRNA splicing. Similarly, SRSF1 coimmunoprecipitated 
with unmodified HDAC2, but not with HDAC2-S394ph 
in MCF7 and HEK293 cells (Supplementary Figure SI). 

As splicing regulator HuR proteins interact with 
HDAC2 and inhibit its activity (17), we measured 
the HDAC activity in the SRSF1 immunoprecipitate. 
Figure IB shows that the HCT116 cell lysate fraction 
immunoprecipitated by antibodies against SRSF1 had 
HDAC activity, which was sensitive to the HDAC inhibi- 
tor TSA or to heat inactivation. Furthermore, the 
addition of a recombinant GST-SRSF1 fusion protein to 
HDAC2 complexes immunoprecipitated from a HCT116 
cell lysate with antibodies against HDAC2 did not affect 
HDAC activity (data not shown). Thus, SRSF1 did not 
inhibit HDAC activity. 

HDAC1 and HDAC2 co-chromatin immunoprecipitate 
with SRSF1 along the body of transcribed genes 

Following transcription initiation, RNAPII pauses until 
it becomes elongation-competent when its CTD is 
phosphorylated at serine 2 (RNAPIIS2ph). SRSF1 is re- 
cruited to RNAPIIS2ph, as soon as the transition to pro- 
ductive elongation occurs (40). In agreement with this 



result, SRSF1 was shown to accumulate along the body 
of the FOS gene only upon induction (41). Similarly, we 
found that SRSF1, HDAC1 and HDAC2 were recruited 
to the body of the FOSL1, MCL1 (Supplementary Figure 
S2) and TFF1 (data not shown) genes upon induction. To 
determine if HDAC1 and HDAC2 are co-chromatin 
immunoprecipitated with SRSF1 to the body of 
transcribed genes, we studied three different genes in 
two cell lines. Serum-starved HCT116 cells were 
stimulated with the phorbol ester, TPA, to activate the 
MAPK pathway and induce the expression of immedi- 
ate-early genes (i.e. FOSL1) (35). To determine the co- 
occupancy of HDAC2/SRSF1 or HDAC1/SRSF1 on 
FOSL1 exons 1 and 4, sequential ChIP assays were per- 
formed, using the dual crosslinking high resolution ChIP 
assay sequentially with antibodies against HDAC2 or 
HDAC1 and SRSF1. Prior to TPA induction, there was 
no HDAC2/SRSF1 or HDAC1/SRSF1 associated with 
these exons, but following TPA induction, HDAC2/ 
SRSF1 or HDAC1/SRSF1 loaded onto the FOSL1 gene 
(Figure 1C and Supplementary Figure S3 A). The expres- 
sion of the MCL1 gene is also increased in response to the 
MAPK pathway activation (42). Accordingly, we 
observed an increased HDAC2 or HDAC1 and SRSF1 
co-occupancy on MCL1 exons 1 and 3, upon TPA induc- 
tion (Figures ID and Supplementary Figure S3A). 
Likewise, induction of the TFF1 gene in MCF7 cells by 
either estrogen or TPA (43,44) resulted in the co-occu- 
pancy of HDAC2 and SRSF1 along the body of 
the gene (Supplementary Figure S3C). Thus, the co-occu- 
pancy of HDAC 1 or HDAC2 and SRSF1 along the body 
of transcribed genes was RNA-dependent, and occurred 
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whether splicing was constitutive (FOSL1 and TFF1) 
or alternative (MCL1). 

SRSF1, HDAC1 and HDAC2 regulate alternative splicing 
of MCL1 RNA 

SRSF1 regulates the alternative splicing of many genes 
(32), including the splicing of MCL1 pre-mRNA, which 
undergoes alternative splicing of exon 2 (45) (Figure 2A). 
Hence, toward the understanding of the mechanisms by 
which HDAC inhibitors alter pre-mRNA splicing, we 
selected MCL1 as a model gene. It has been previously 
shown that a decline in availability of SRSF1 favored the 
production of the pro-apoptotic MCL1 short form 
(MCL1S) transcript (exclusion of exon 2) (45). Similarly 
we found that knockdown of SRSF1 affected the splicing 
of MCL1 transcripts in HCT116 cells (Figure 2B). We 
observed a greater exclusion of alternative exon 2 in 



HCT116 cycling cells transiently transfected with SRSF1 
siRNA compared to cells transfected with scramble 
(nontargeting) siRNA. Four independent sets of data 
showed that the increase of the percentage of MCL1 
exon 2 exclusion (MCL1S transcripts) in SRSF1 
knockdown cells had a strong statistical significance 
(Figure 2B). A value of 84.5 ± 5.3% SRSF1 knockdown 
was determined by immunoblot analysis (Supplementary 
Figure S4A). 

Next, we determined the effect of HDAC inhibition on 
the alternative splicing of MCL1. To control for any off- 
target effect, we used two HDAC inhibitors with different 
selectivities: TSA or apicidin. TSA is considered a pan- 
inhibitor; meaning it inhibits Class I and lib HDACs. 
On the other hand, apicidin is a Class I HDAC inhibitor, 
inhibiting specifically HDAC1, 2 and 3 (46). To minimize 
secondary effects of HDAC inhibitors treatment, such as 
cell-cycle arrest or apoptosis, we limited treatment times 
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to a maximum of 2 h. In cycling HCT1 16 cells treated with 
TSA or apicidin for only 2 h, a shift in the splicing of the 
MCL1 pre-mRNAs was observed in favor of MCL1S 
(Supplementary Figure S4B and Figure 2C). We 
repeated our studies with a natural HDAC inhibitor, 
butyrate at a physiologically relevant concentration 
(47,48) (Supplementary Figure S4C). With the use of 
three different HDAC inhibitors which are structurally 
unrelated and have different spectrum of substrate 
specificities, similar results were obtained (Figure 2C, 
Supplementary Figures S4B and S4C). Moreover, we 
tested the expression level of SRSF1 following treatment 
with HDAC inhibitors. Immunoblot analysis showed that 
the change in splicing resulting from HDAC inhibition 
was not due to a reduced expression of SRSF1. Up to 
24 h treatment with TSA or apicidin did not affect 
SRSF1 protein levels (Supplementary Figure S4D). As 
SRSF1 is modified by acetylation (49), we probed its 



acetylation level with or without 2 h TSA or apicidin treat- 
ment. Reciprocal coimmunoprecipitation assays showed 
that the acetylation level of SRSF1 was not affected by 
HDAC inhibitors, while the histone H3 acetylation level 
was markedly increased (Supplementary Figure S4E). 
Thus, we were able to rule out an indirect effect of TSA 
or apicidin on MCL1 splicing through differences in 
SRSF1 levels or SRSF1 acetylation. As a complementary 
approach and to specifically study the outcome of loss of 
HDAC1 and/or HDAC2 activity, we determined the 
effects of HDAC1, HDAC2 or HDAC1/HDAC2 
knockdown on the MCL1 splicing pattern. A representa- 
tive immunoblot analysis of HDAC1 and/or HDAC2 
knockdown in HCT116 cells is shown in Supplementary 
Figure S4F. HDAC1 knockdown was determined to be 
78 ± 1% and 76 ± 2% in HDAC1 and HDAC1/2 
knockdown HCT116 cells, respectively, while HDAC2 
knockdown were 84 ± 10% and 85 ± 12% in HDAC2 
and HDAC 1/2 knockdown cells, respectively. We show 
that knockdown of either HDAC1 or HDAC2 or both 
led to a strongly significant change in splicing in favor 
of exon 2 skipping (Figures 2D and Supplementary 
Figure S4G). The immunoblot in Supplementary Figure 
S4F shows that the SRSF1 level was not affected by 
HDAC knockdown. These results indicate that HDAC1 
and HDAC2 play a role in splice site selection by a mech- 
anism other than altering expression or acetylation levels 
of SRSF1. 

HDAC inhibitors affect the MCL1 gene occupancy 
by RNAPII 

The current literature supports the kinetic model where 
variation in RNAPII elongation rate or processivity 
modulates the splice site recognition efficiency. For 
genes with alternative exons that have weak splice sites, 
a slowly moving RNAPII results in the inclusion of the 
exons in the RNA (50). To investigate the effect of HDAC 
inhibitors (apicidin, TSA, butyrate) on the MCL1 gene 
occupancy by the elongating form of RNAPII 
(RNAPIIS2ph), HDAC2, SRSF1 and acetylated histone, 
we performed dual crosslinking high resolution ChIP 
assays. Following a 2 h treatment of cycling HCT116 
cells with the HDAC inhibitors, there was a reduced as- 
sociation with RNAPIIS2ph, HDAC2 and SRSF1 with 
the exon 2 nucleosome relative to nucleosomes at exon 1 
and 3 (Supplementary Figures S5, S7 and S8). Also there 
was an increase in H3 and H4 acetylation of the exon 2 
nucleosome (Supplementary Figures S6, S7 and S8). 
Meanwhile, results of H3 ChIP assays showing a 
uniform association of H3 along the MCL1 gene with 
or without apicidin (Supplementary Figure S5) are con- 
sistent with ChIP assays being performed on mono- 
nucleosomes obtained by MNase digestion of chromatin 
(Supplementary Figure S5 and S7). 

As dynamic acetylation results from the balance of 
opposing activities of HDACs and KATs, we turned our 
attention to KATs. A genome-wide mapping of several 
KATs revealed that p300 and CBP (KAT3A and B) as 
well as TAF1 (KAT4) were associated with promoter 
regions, while KAT2B (PCAF), KAT5 (TIP60) and 
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KAT8 (MOF) were located with transcribed regions (15). 
Moreover, KAT2B and KAT7 (HBOl) have been 
implicated in the elongation phase of transcription 
(51-53). Both KAT2B and KAT7 are H3 and H4 modify- 
ing enzymes (54). Hence, their distribution along the 
MCL1 gene body was studied. The distributions of 
KAT2B and KAT7 were opposed to that of HDAC2, 
with an increased loading on exon 2 (Supplementary 
Figure S5 and S7). 

Next we determined whether SRSF1 had a role in 
modulating the chromatin structure and/or the RNAP11 
elongation rate, which could lead to the change in splicing 
of MCL1 transcripts. To this end, we compared the dis- 
tributions of RNAP11, HDAC2 and acetylated histones 
along the MCL1 gene in scramble and SRSF1 knockdown 
cells by ChIP assays. Under these conditions, no obvious 
changes in HDAC2 association, RNAPII distribution or 
acetylated histones were detected (Supplementary Figure 
S9). These results suggest that SRSF1 is involved in alter- 
native splicing regulation of MCL1 in a mechanism inde- 
pendent of the modulation of histone acetylation. 

Effect of HDAC inhibitors on MCL1 splicing is amplified 
upon TPA stimulation of serum-starved cells 

Considering the observed consequences of HDAC inhib- 
ition on chromatin structure surrounding exon 2, the 
effects on splicing (Figure 2C) were rather small. This dis- 
crepancy was also observed previously and was interpreted 
as interference from the mRNA synthesized and spliced 
before HDAC inhibition (16). To alleviate this problem, 
the authors treated the cells with DRB, a reversible inhibitor 
of RNAPII, prior to HDAC inhibition. In this study, we 
exploited the transcriptional response of the MCL1 gene 
downstream of the MAPK pathway (42) to study the 
effects of HDAC inhibitors on splicing during a 
synchronized induction of MCL1 expression. Serum- 
starved HCT116 cells, pretreated or not for 30min with 
apicidin or TSA, were stimulated with TPA. TPA stimula- 
tion of MCL1 expression in control cells resulted in slight 
increase in exon 2 exclusion [see apicidin (— ), TPA 60min in 
Figure 3A; TSA (— ), TPA 30 and 60 min in Supplementary 
Figure SI OA]. With an increase in transcription initiation at 
the MCL1 promoter, a higher density of elongating 
RNAPII was expected along the body of the gene. 
Indeed, RNAPIIS2ph ChIP assays revealed an increase in 
the number of RNAPII molecules along the body of the 
MCL1 gene [see apicidin (— ), TPA 60 min in Figure 3C; 
TSA (— ), TPA 30 and 60 min in Supplementary Figure 
SI 0B]. This increased occupancy along the MCL1 gene 
was also observed for SRSF1, HDAC2, KAT2B and 
KAT7 (Figure 3C and Supplementary Figure S10B), sug- 
gesting that all three proteins were associated with the 
elongating RNAPII complex. MCL1 TPA stimulation was 
also accompanied by a global H3 and H4 acetylation (as 
well as H3K9ac, H3K14ac, H4K5ac and H4K8ac) increase 
along the body of the MCL1 gene (Figures 4A, 
Supplementary Figures Sll and SI 2). 

Besides interfering with dynamic acetylation of 
histones, HDAC inhibitors upregulate MCL1 transcrip- 
tion (55). Thus, we analyzed the MCL1 splicing pattern 



in serum-starved HCT116 cells submitted to a 30 min 
preincubation with apicidin or TSA, just prior to TPA 
induction. Similar to the TPA-induced control cells, 
apicidin- or TSA-treated cells had a slightly greater level 
of MCL1 exon 2 exclusion than control cells [see apicidin 
(+), TPA Omin in Figure 3A; TSA (+), TPA Omin in 
Supplementary Figure SI OA]. Apicidin- or TSA-treated 
cells also exhibited an increased occupancy of 
RNAPIIS2ph, SRSF1 and HDAC2 along the body of 
the MCL1 gene compared to untreated cells, reflecting 
the HDAC inhibitor-induced upregulation of transcrip- 
tion initiation. However, at the same time, KAT2B and 
KAT7 started to associate specifically with exon 2. 
Moreover, a parallel exon 2-specific H3 and H4 acetyl- 
ation increase was observed (Figures 4A, Supplementary 
Figures Sll and SI 2). So, after 30 min of inhibition of 
HDAC activity in serum-starved HCT116 cells, KAT2B 
and KAT7 have started to target the nucleosome pos- 
itioned on the alternative exon. 

Following TPA stimulation, apicidin and TSA had a 
profound impact on splicing, resulting in a marked 
increase of MCL1 exon 2 exclusion [see apicidin (+), 
TPA 60 min in Figure 3A; TSA (+), TPA 30 and 60 min 
in Supplementary Figure SI OA]. In the meantime, while 
RNAPIIS2ph on exon 1 increased as a consequence of 
TPA induction, the occupancy over exon 2 dropped dra- 
matically, suggesting a sharp acceleration of the elong- 
ation process over exon 2. At exon 3, RNAPIIS2ph 
occupancy was similar to those on exon 1 [see apicidin 
(+), TPA 60 min in Figure 3C; TSA (+), TPA 30 and 
60 min in Supplementary Figure SI 0B]. Moreover, the dis- 
tribution of SRSF1 and HDAC2 on the body of the 
MCL1 gene was consistent with these proteins being 
associated with RNAPIIS2ph. On the other hand, 
KAT2B and KAT7 levels rose sharply at exon 2 
(Figures 3C and Supplementary Figure SI 0B), as did 
overall H3 and H4 acetylation, as well as H3K9ac, 
H3K14ac, H4K5ac, H4K8ac (Figures 4A, 
Supplementary Figures Sll and S12). It should be noted 
that changes at exon 2 were gradual changes occurring 
over time. This is particularly noticeable in the TSA-in- 
hibition experiment (Supplementary Figures SI 0B and 
SI 2). As a control, the total H3 distribution along the 
body of the MCL1 gene was analyzed. It was unaffected 
by TPA induction with or without apicidin or TSA treat- 
ment (Figures 3C and Supplementary Figure S10B). 

Previous studies have demonstrated that H3 
trimefhylated at K4 was preferentially engaged in 
dynamic acetylation (15,56). In ChlP-sequencing studies, 
which determined the genomic distribution of H3K4me3 
in HCT1 16 cells, we observed intense H3K4me3 at exon 2 
of the MCL1 gene (Figure 4B). This observation is in ac- 
cordance with the selective histone acetylation of the exon 
2 nucleosome, suggesting that the H3K4 methylated state 
of the exon 2 nucleosome has a role in recruiting KATs. 

To recapitulate the above results, the TPA induction of 
MCL1 expression was accompanied by an overall increase 
in density of RNAPIIS2ph and a slight increase in alter- 
native exon 2 exclusion. On the other hand, the TPA in- 
duction of MCL1 expression with HDAC inhibition led to 
gradual important local changes on exon 2 chromatin, 
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Figure 3. HDAC inhibition favors exclusion of alternative exon 2 upon TPA induction of MCL1 gene. (A) MCL1 alternative exon 2 splicing was 
analyzed in TPA-stimulated serum-starved HCT116 cells treated or not with 150nM apicidin prior to TPA induction. (B) Schematic representation of 
the amplicons generated by real-time PCR analyses of ChIP DNA. (C) ChIP assays were performed on DSP- and formaldehyde-crosslinked 
mononucleosomes prepared from TPA-stimulated serum-starved HCT116 cells treated or not with 150nM apicidin for 30 min prior to TPA induction. 



resulting in exclusion of exon 2 from the pre-mRNA. 
Increased acetylation of H3 and H4 on nucleosome 
residing on exon 2 and targeting of KAT2B to exon 2 
preceded the apicidin- or TSA-induced RNAPIIS2ph 
reduced occupancy over exon 2 [see apicidin (+), TPA 
Omin in Figures 3C and 4A; TSA (+), TPA Omin in 
Supplementary Figures SI OB and Sll]. These results 
suggest that the destabilization of the nucleosome over 
exon 2 leads to an increased elongation rate of this exon. 

HDAC1, HDAC2 and KAT2B are associated with both 
the pre-mRNA and chromatin 

ChIP assays have shown that the association of several 
splicing factors (SR proteins and HuR) to transcribed 



genes was sensitive to RNase treatment, indicating that 
these proteins bind more strongly to the pre-mRNA 
than to elongating RNAPII (17,41). To determine 
whether this was the case for SRSF1, we performed 
ChIP assays on DSP- and formaldehyde-crosslinked 
mononucleosomes prepared from lysates of TPA- 
stimulated serum-starved HCT116 cells, which were 
treated or not with RNase A prior to the immunopre- 
cipitation step. Our results show that RNase treatment 
of dual crosslinked mononucleosomes reduced the associ- 
ation of SRSF1 with the body of the MCL1 gene 
in serum-starved and TPA-stimulated HCT116 cells 
(Figure 5A), suggesting that SRSF1 was bound to the 
pre-mRNA. As HDAC2 is corecruited with SRSF1 to 
the body of transcribed genes, we repeated our ChIP 



Nucleic Acids Research, 2014, Vol. 42, No. 3 1665 



g 7 



CD 0 
?1 



c 7 



S3 
11 



B 



TPA 0 TPA 60 

H3K14acChlP 




Apicidin (-) Apicidin (+) 

H4K5acChlP 



c 7 
cu 

15 
o 

cu 3 
■a 
o 1 



H3K9ac ChIP 













Apicidin (-) 


Apicidin (+) 



H4K8ac ChIP 




Apicidin (-) Apicidin (+) 



Apicidin (-) 



0\ <c> 0> 0 

Apicidin (+) 



13 



I 7 



MCL1 











li 


I 


L 1 


111 I 1 1 1. 





150551.8 


150549.9 


150548.1 








E1 


E2 


E3 



KB 



Figure 4. Specific enrichment of H3 and H4 acetylation and methylation over MCL1 exon 2. (A) HDAC inhibition specifically increases H3 and H4 
acetylation over exon 2 upon TPA induction of MCL1 gene. ChIP experiments were performed as described in Figure 3. (B) ChlP-sequencing determined 
H3K4me3 signals along the MCL1 gene were displayed by Partek Genomic Suite. The Y-axis is the signal intensity of mapped sequence reads. 



assay with antibodies against HDAC2. Figure 5A shows 
that HDAC2 association with MCL1 was markedly 
reduced upon RNase digestion, suggesting that HDAC2 
interacted with the MCL1 pre-mRNA. 

It was shown that SRSF1 interaction with 
RNAPIIS2ph was also mediated by the pre-mRNA 
given that it was sensitive to RNase as opposed to a 
direct protein-protein interaction, which would be 
RNase-insensitive. We set out to confirm this finding in 
HCT116 cells and find out if HDAC2 also formed 
a RNase-sensitive cotranscriptional complex with 
RNAPIIS2ph. Thus, we carried out reciprocal coimmuno- 
precipitation assays using HCT116 cell lysates treated or 
not with RNase A. Figure 5B shows that SRSF1 inter- 
action with RNAPIIS2ph was RNA-dependent as SRSF1 
did not coimmunoprecipitate with RNAPIIS2ph in the 
presence of RNase A, and vice versa (Supplementary 
Figure SI 3). Similar results were obtained for HDAC2 



and RNAPIIS2ph interactions (Figures 5B and 
Supplementary Figure SI 3). 

To confirm the association of SRSF1 and HDAC2 with 
the MCL1 mRNA, we used the CLIP (ultraviolet 
crosslinking and immunoprecipitation) method to 
analyze their in situ direct binding to RNA. We also 
investigated the interactions of HDAC1 and KAT2B 
with RNA, as KAT2B had been reported to associate 
with RNAPII, in complex with hnRNP U and actin 
(52). RNP complexes in cycling HCT116 cells, treated or 
not for 2h with apicidin or TSA, were covalently 
crosslinked by UV exposure only or by DSP and UV 
exposure, and RNA complexes were immune-purified 
with anti-SRSFl, anti-HDACl, anti-HDAC2 or anti- 
KAT2B antibodies. In UV-treated cells, RT-PCR 
analysis showed a significant enrichment of MCL1 exon 
1, exon 2 and exon 3 mRNAs in SRSF1 immunopre- 
cipitation versus an irrelevant IgG control, but the 
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enrichments in HDAC1, HDAC2 or KAT2B immunopre- 
cipitations were minimal for all three exon mRNAs 
(Figure 6). On the other hand, when RNP complexes 
were crosslinked by DSP and UV exposure, there was a 
significant enrichment of MCL1 exon 1 , exon 2 and exon 3 
mRNAs in each of the immunoprecipitations versus an 
IgG control. The results were independent of HDAC 
activity, as they were not affected by the treatment of 
cells with HDAC inhibitors apicidin or TSA (Figure 6). 
The detection of significant HDAC1, HDAC2 and 
KAT2B association with the MCL1 mRNA being depend- 
ent on dual crosslinking indicates that the interactions of 
these proteins with the mRNA are indirect and are 
mediated by splicing factors. 

Together, these results are consistent with the idea that 
HDAC1, HDAC2 and KAT2B are associated with the 
pre-mRNA while catalyzing the dynamic acetylation of 
MCL1 chromatin. 



DISCUSSION 

Nucleosomes are nonuniformly distributed along the body 
of transcribed genes and are enriched on exons relative to 
introns. Notably, the average length of human internal 
exons is very close to the length of DNA within a nucleo- 
some (147 bp). This marked correspondence also applies 
to six other tested metazoans (57,58), suggesting a 
conserved role for the nucleosome in exon definition and 
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Immunoprecipitations with indicated antibodies were performed on 
UV light-exposed or dual DSP and UV light crosslinked RNP 
complexes isolated from HCT116 cycling cells. Reverse transcription 
PCR measurements were normalized to the value obtained with 
IgG control antibodies. 
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splicing regulation. Our studies demonstrate that the nu- 
cleosome positioned over the alternative exon 2 of the 
MCL1 gene is highly dynamic with regard to histone 
acetylation, a state which greatly impacts on splicing 
decision. 

Our mass spectrometry data showed that HDAC1 and 
HDAC2, more accurately nonphosphorylated HDAC2, 
coimmunoprecipitated with numerous splicing factors. 
Interestingly mass spectrometry analyses of the 
spliceosome had identified HDAC2 as a component of 
this complex (59). SR proteins associate with chromatin 
in a transcription-dependent manner (41), and we show 
that HDAC1 and HDAC2 are co-chroma tin immunopre- 
cipitated with SRSF1 along the body of genes upon tran- 
scription induction, regardless of whether the transcripts 
of these genes are constitutively or alternatively spliced. 
The presence of SRSF1 in the HDAC1 and HDAC2 
immunoprecipitate did not inhibit HDAC activity as the 
SRSF1/HDAC1 and SRSF1/HDAC2 immunoprecipitates 
had deacetylase activity. Of note, we repeated the in vitro 
experiment showing that HuR proteins inhibit the HDAC 
activity in a HDAC2 immunoprecipitate from HeLa 
nuclear extracts (17), but contrary to what was published, 
we found no impact of the HuR protein on HDAC 
activity. It was proposed that SRSF1 is recruited to 
transcribed genes by RNAPII upon stimulation of elong- 
ation, which occurs with the LARP7 (P-TEFb)-mediated 
phosphorylation of the CTD of RNAPII at Ser 2 (40). It is 
most likely that HDAC1 and HDAC2 are recruited at the 
same time, in association with splicing factors. This model 
is supported by our data showing that HDAC2 inter- 
actions with numerous splicing factors are unaffected by 
RNase treatment. Furthermore, SRSF1 knockdown did 
not affect the recruitment of HDAC1 and HDAC2 to 
the body of the MCL1 gene or the TFF1 gene upon its 
induction (data not shown). Also the increased MCL1 
exon 2 exclusion resulting from a SRSF1 knockdown 
did not affect MCL1 exon 2 acetylation or the RNAPII 
distribution along the MCL1 gene (60), suggesting that 
SRSF1 affects MCL1 splicing through a mechanism 
which is not dependent on modulation of histone acetyl- 
ation. Since HDAC1 and HDAC2 are associated with so 
many splicing proteins, it is not surprising that their re- 
cruitment to the body of a transcribed gene is not affected 
by the knockdown of one of them. Vice-versa, HDAC1 
and/or HDAC2 knockdown did not affect SRSF1 recruit- 
ment to the TFF1 gene. 

SR proteins interact with the nascent pre-mRNA, 
RNAPIIS2ph and chromatin. SR proteins interactions 
with RNAPIIS2ph are RNase-sensitive, indicating that 
they are RNA mediated. On the other hand, SR 
proteins interactions with chromatin were reduced by 
RNase treatment, but not abolished, suggesting that 
these interactions are partially mediated by RNA (41). 
Indeed, in the case of SRSF1, interaction with chromatin 
also occurs via the tail of nucleosomal histone H3 (61). In 
this study, we show for the first time that HDAC1 and 
HDAC2 also reside on the pre-mRNA. In fact, the asso- 
ciation of HDAC2 with the body of the MCL1 gene was 
particularly sensitive to RNase, more so than the associ- 
ation of SRSF1. Likewise, HDAC2 association with 



RNAPIIS2ph was RNA-dependent, indicating that 
HDAC2 was not bound directly to RNAPIIS2ph. Our 
results demonstrating that HDAC1 and HDAC2 are 
associated with the emerging RNA transcripts are sup- 
ported by our previous ChIP experiments, in which we 
needed to use dual crosslinking to monitor HDAC1 and 
HDAC2 along the body of genes; formaldehyde alone 
worked poorly (12). Likewise, a genome wide-mapping 
study of KATs and HDACs applied dual crosslinking 
(15). In agreement with previous results revealing that 
KAT2B is associated with RNAPII, in complex with 
hnRNP U and actin (52), we show that KAT2B is 
bound to the pre-mRNA. We propose that HDAC1 and 
HDAC2, in concert with KAT2B, and other KATs, 
catalyze nucleosomal dynamic acetylation by acting at 
the RNA level. 

Evidence is accumulating to show the importance of 
nuclear RNA in maintaining the structure of transcribed 
chromatin as well as the nuclear location of transcribed 
genes (62-65). It is quite possible that nuclear RNA 
associated with regulatory and coding regions of 
transcribed and silent chromatin domains serves as a 
platform for several chromatin-modifying enzymes. 

It is interesting to note that the ENCODE ChlP- 
sequencing data show that H3K4me3 peaks at the 
MCL1 exon 2 in many cell types (66,67). There are 
several mechanisms by which the exon 2 nucleosome 
H3K4me3 mark recruits KAT2B and KAT7. KAT2B is 
a component of the SAGA and ATAC complexes which 
contain the SAGA-associated factor 29, which binds to 
H3K4me3 (68,69). Similarly the ING4 subunit of the 
KAT7 complex binds to H3K4me3 (70). The histone 
acetylation activity of the KATs is countered by the 
HDACs. It should be noted that in cycling cells, the 
steady state levels of H3K9ac and H3K14ac at exon 2 
are decreased compared with the rest of the MCL1 gene 
body. This observation suggests that HDAC activity is 
locally increased relative to KAT activity on exon 2. 

Our data show that the recruitment of KATs to exon 2 
and increased acetylation of exon 2 nucleosome precedes 
locally decreased RNAPII occupancy, implying that 
increased nucleosomal acetylation contributes to faster 
elongation. This conclusion is supported by an in vitro 
transcriptional elongation experiments showing that 
acetylation of H3 and H4 tails overcomes the inhibitory 
effect of nucleosomes (71). In situ experiments have also 
demonstrated that histone acetylation is required to 
maintain the unfolded structure of the transcribed nucleo- 
some (72). It is assumed that histone acetylation facilitates 
subsequent rounds of transcription elongation (73). 

Upon HDAC inhibition or knockdown (HDAC1 and/ 
or 2), exon 2 nucleosome acetylation markedly increases 
due to the H3K4me3 localized KAT activity. Within 
30 min of HDAC inhibition, an enhanced recruitment of 
KAT2B and KAT7 to MCL1 exon 2 was observed. 
Recruitment of the KAT may escalate through their 
binding to acetylated H3 and H4 through the enzyme's 
bromodomain (74). By binding to the product of its 
activity, KAT2B occupancy could be self-reinforcing. 

SR proteins regulate alternative splicing by binding to 
exonic splicing enhancers and enhancing Ul snRNP and 
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U2 snRNP recruitment, thus promoting recognition of 
exons with suboptimal splice sites (41). In agreement 
with this mode of action, we and others (45) found that 
SRSF1 knockdown resulted in increased MCL1 exon 2 
exclusion. Enzymatic inhibition of HDACs (three 
HDAC inhibitors, TSA, apicidin and butyrate that are 
structurally unrelated and have different substrate 
spectra) or siRNA-mediated knockdown of HDAC1 
and/or HDAC2 also resulted in increased MCL1 exon 2 
exclusion, without affecting SRSF1 levels. Of these 
HDAC inhibitors, butyrate produced by fermentation of 
dietary fiber by colonic bacteria could reach sufficient 
levels (0.5-5 mM) to alter MCL1 splicing in colonic cells 
(47,48). Following HDAC inhibition, the hyperacetylation 
of MCL1 exon 2 nucleosome and increased elongation 
rate may result in an altered localization and decreased 
recruitment of SRSF1 and other splicing factors to the 
MCL1 transcript (23). These events could explain 
the loss of association of SRSF1 and HDAC 1/2 with the 
hyperacetylated MCL1 exon 2 nucleosome in our ChIP 
assays (23). Thus, although acting through different mech- 
anisms, HDAC inhibitors and SRSF1 knockdown prevent 
the SRSFl-mediated exon 2 inclusion (23). 

Future studies will determine the role of H3K4 methy- 
lation and dynamic histone acetylation in regulating alter- 
native RNA splicing. We propose that HDAC inhibitors 
will alter pre-mRNA splicing when nucleosomes placed at 
regions containing alternative exons are modified by 
H3K4 methylation. The nucleosomes marked by H3K4 
methylation are preferentially engaged in dynamic acetyl- 
ation (15,56). Hence HDAC inhibition may increase 
elongation rates at the alternative exon and alter the as- 
sociation of splicing factors to the transcript, resulting in 
exon exclusion or inclusion (23). To conclude, our results 
highlight the roles of HDACs and KATs in regulating the 
interface between chromatin organization and alternative 
splicing. 
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